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ABSTRACT 



Solutions have been obtained for the temperature and moisture distribu- 
tions within the atmospheric boundary layer of an axisymmetric hurricane 
model. The intensity of the hurricane is related to the equivalent potential 
temperature gradient produced by a balance of heat sources and advection 
within the boundary layer. Solutions are obtained using the bulk aerodynamic 
transport equations or applying a two-layer, baroclinic boundary layer model 
by Cardone (1969). Equilibrium maximum wind speeds vary from minimal hurri- 
cane force to about 80 m/sec for fixed ocean temperatures between 27C and 31C. 
Observations of the oceanic heat loss in actual storms are necessary to estab- 
lish the ratio of heat flux to momentum transfer in the high wind speeds 
treated by the model. Various applications of the model are proposed, with 
an ultimate goal of a time-dependent simulation of hurricane-ocean coupling. 

INTRODUCTION 

A number of axisymmetric vortex models (e.g., Yamasaki, 1968a and 1968b; 
Ooyama, 1969; Rosenthal, 1971) have been developed to describe the evolution 
of the intense wind region characteristic of hurricanes. The specified 
initial conditions are mainly designed to insure low-level inflow, so that 
the vortex will eventually become unstable, and are not related to observa- 
tions. As a result, a comparatively long "organizational period" is required 
before deepening of the vortex ensues. A primary result of these models is 
the very realistic simulation of the wind structure near the center of the 
hurricane. Because the initial conditions are arbitrarily specified, rather 
than from real data, the intermediate stages of the vortex simulation may not 
be realistic. These models do not allow forcing by external circulations. 



Instead, the philosophy is to displace the exterior boundaries of the model 
as far as necessary to assure that the boundary conditions will not affect 
the solution near the center. A complete assessment of external circula- 
tions, and other asymmetric effects, on hurricane simulation must await fur- 
ther development (see Anthes et al, 1971) of three-dimensional models started 
with real initial data. 

The objective of this paper is to describe a simplified model of the 
properties within the hurricane boundary layer. Although the ultimate pur- 
pose is to simulate the time variation of these properties due to changes in 
sea-surface temperatures, in the present paper the emphasis will be on the 
adjustment to different fixed sea-surface temperatures. It is well known 
that tropical storms develop into hurricanes only over warm tropical oceans. 
Indeed the oceanic surface layers must provide the heat necessary to produce 
the unique warm core of the mature hurricane. The preferential release of 
the heat in deep convective towers near the hurricane eye results in central 
temperatures much higher than the environment. To a first approximation, 
the heat gain within the boundary layer may be linearly related to the 

tangential wind speed at the top of the layer (vg), but the frictional loss 

2 

of momentum to the ocean increases as Vg . Thus a given radial profile of 
wind speed essentially specifies the heat gain which sustains the hurricane, 
as well as the frictional loss which tends to offset further increases in 
wind speed. On this basis, Riehl (1963) developed relations for the unique 
wind profile (and associated thermal structure) for which these offsetting 
tendencies are in balance; that is, for the conditions in which a steady 
state hurricane should result. However, an actual storm may not reach its 
full potential, or may be in steady state for only a short time, because the 
characteristics of the ocean and surrounding environment vary in space and 
time. A prime motivation of this study was to develop a model capable of 
simulating the symmetric component of actual storms. 
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DESCRIPTION OF THE MODEL 



The empiricism of the model is based on three assumptions. First, the 
wind profile at the top of the boundary layer is specified by the relation 

1/2 

Vgr = constant , (1) 

outside the radius of maximum wind speed (r^). Between the vortex center and 
r^ the winds are assumed to obey Vgr ^ = constant, which represents solid- 
body rotation. Justification of (1) was given by Gray and Shea (1973) who 
found an exponent of 0.47 + 0.30 in Vg profiles from hurricanes in various 
stages of development. Rather than using individual Vg profiles, Riehl (1963) 
averaged profiles around the storm and found the exponent (1/2) applied in 
many cases, particularly in quasi-steady state hurricanes. Some of the addi- 
tional variation which Gray and Shea found must be due to transient features 
in the wind field as the storm adjusts to varying internal or external condi- 
tions. In simulating the free-atmosphere as a series of ’’equilibrium states” 
represented by (1), these transient features are assumed to be of less 
importance than the symmetrical circulation. It should be emphasized that 
the coefficient in (1) was determined from aircraft data only within about 
150 km of the center. This is the region in which the symmetry assumption is 
most justifiable, and, fortunately, is the primary region of interest in 
evaluating the effects of the hurricane. Lack of systematic observations at 
larger radii prevent an accurate description of the outer region in which 
the external effects are likely to be manifested in actual hurricanes. There- 
fore, beyond r^ [defined in (5)] the exponent in (1) was adjusted to always 
reduce the wind speed at 2r^ to 5 m sec Specifying a different exponent 

in the outer region had an insignificant effect on the properties near the 
center compared to experiments in which (1) was applied throughout the region. 
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The second empirical formula relates the maximum wind speed to the in- 
crease in the equivalent potential temperature of the boundary layer (‘^6^) 
due to the oceanic heat source. Such a relation was proposed by Riehl (1963), 
through the use of the cyclostrophic assumption and observations relating 
the surface pressure at the radius of maximum wind to 69^, 

6p = - 2,5 66^ , (2) 

where the 6 symbols indicate deviations from p = 1005 and 6^ = 350K. In the 
present model, a gradient wind expression 
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was used instead of the cyclostrophic assumption. Although Ooyama (1969) 
showed that the inertial term is important near the surface, this term may be 
dropped if (3) is averaged over the depth of the boundary. Upon integration 
of (3) between the radii over which (2) applies, the tangential component 
(Vq ) at the radius of maximum wind (r.) is 
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The dominant term in this expression is that involving the difference in 0^ 
at r. and at some radius r where p = 1005 and 0 = 350K, which follows 

from (2). Bell and Tsui (1973) found a similar relation between pressure 
and 6^ changes for tropical cyclones near the Asian coast, but with a co- 
efficient of 2,25. Because of the difficulty in establishing the desired 
surface pressure at r^, and only a 10^ percent difference, the coefficient 
in (2) was adopted. It is clear from (4) that the radius of maximum wind 
(r^) should coincide with the maximum 0^ within the boundary layer. In the 
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present model, was chosen as the radius at which 0^, averaged over five 
grid increments (Ar = 3km), was a maximum. This assumption makes r^ coin- 
cident with the eyewall clouds, because the boundary layer air parcels 
with the largest heat and moisture content would be expected to be most 
favorable for deep convection. Shea and Gray (1973) show that r^ tends to 

be coincident with the eyewall in the most intense hurricanes, but r. is 

1 

larger than the eyewall radius in weaker storms. The radius (r^) at which 
p = 1005 and 0^ = 350K may be taken to represent an external influence, and 
must be defined in terms of other parameters. As the p and 0^ gradients are 
small outside the central region, the inaccuracy in specifying r^ is a 
second-order effect compared with r^. For convenience we choose 



r 

o 



,2 , 1/2 
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(5) 



1 

which is approximately the radius at which Vq tends to zero in the outflow 
layer of the storm, if the air originating in the maximum wind region con- 
serves absolute angular momentum. Since r^ is about 200km for f^(10 ^)sec 

and typical values of Vg and r^, this again represents the region of appli- 

i 

cation of the model. A schematic representation of the wind structure in 

the atmospheric vortex is given in Figure 1. 

The above relations represent an empirical model of the wind speeds in 

hurricanes as a function of the 0 distribution within the boundary layer. 

e 

Rather than determine the maximum 0^ from observations as proposed by Riehl 
(1963), the 0^ fields will be calculated from time -dependent equations for 
the potential temperature (0) and specific humidity (q). 



9^ = 6 exp 

P 

jg _ S(^„e) f5s_ 
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( 6 ) 



5 



( 7 ) 



5q _ ^ (rv^q ) ^ (^q ) 

dt rdr Sp LAp 



In these equations Cp and L are the specific heat and latent heat of vapori- 
zation. Only the symmetrical circulation is considered, so that the radial 
(v^) and vertical (cd) components of motion are related by 






Srv 



rSr ’ 



( 8 ) 



with the further approximation that oo = 0 at the surface. Radiative exchanges 
are assumed to be negligible and the sensible (Q^) and latent (Q^) heat 
fluxes are assumed to decrease to zero within Ap = 100 mb of the surface. 

The actual vertical heat fluxes through the top of the boundary layer are 
determined by the mean flow represented by o) and convective towers. Normally 
the convective heat flux is parameterized in terms of the mean flow; however, 
this resulted in supersaturation in the present model, as was the experience 
of Rosenthal (1971). The third important empirical aspect of the model was 
the formulation of a critical relative humidity profile, such that the 
moisture in excess of the critical value was arbitrarily removed from the 
boundary layer. The convective heat flux O30 was neglected o 

Because the radial gradient of 0^ is essentially determined by the 
moisture gradient, the specification of the critical relative humidity profile 
was a primary factor in determining the intensity of the model hurricane. 

As a first attempt, the critical relative humidity for the air in the exterior 
region was specified as the value which would produce a lifting condensation 
level (LCL) of 60 mb above the surface. A decrease in the typical cloud base 
height toward the center is commonly observed in hurricanes, although the 
magnitude of the decrease has not been documented as a function of maximum 
wind speed. In this case, a linear decrease to a LCL of 30 mb at the eyewall 
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was assumed. Profiles other than a linear decrease were also used, but 
these did not change the results. These LCL values corresponded to external 
and internal critical relative humidity values of about 77 and 88/S, respec- 
tively, for the range of sea surface temperatures used here. These values 
seem realistic as one might expect the external air to be somewhat drier 
because of subsidence, and the interior air to be more moist than normal 
(about 80^ in the subcloud layer). In the second approach, the critical 
relative humidity values were arbitrarily specified as a linear profile be- 
tween 95^ inside r^ and 75/2 outside r^. These values might be taken as 
representative of the maximum conditions one might expect, and the results 
for this specification will be referred to as the maximum hurricane. A 
95/2 value implies a cloud base of about 150 m near the center. 

BOUNDARY LAYER MODELS 



The surface fluxes Q^, and Tg^ required to close the solutions are 
determined by the boundary layer model. For fixed oceanic conditions, the 
model will tend to steady state as the radial and vertical advection of 
heat and moisture within the boundary layer balance the oceanic sources 
Qg and Q^, respectively. In this formulation the air-sea temperature and 
humidity differences will vary with radius as determined from (6) and (7). 
The first approach to the boundary layer was to assume the bulk aerodynamic 
equations for momentum, heat and moisture fluxes: 



Tn = p V Vq 
0s D s 6s 



(9) 

( 10 ) 
( 11 ) 
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where is the heat /moisture exchange coefficient, is the surface shear- 

ing stress, is the specific humidity corresponding to the saturation vapor 
pressure at the sea surface temperature T^, and the total wind speed (V^) is 
0.7 times the wind speed at the top of the boundary layer. If the shearing 
stress vanishes at the top of the boundary layer, the steady-state, tangential 
equation of motion may be used to determine the mean radial flow within the 
layer. 



V 

r 



SVn 




( 12 ) 



The set (1) - (12) is thus closed for a given sea surface temperature, and 
may be solved as an initial value problem. 

As an alternative to the bulk aerod 3 mamic equation approach, a two-layer, 
baroclinic model developed by Cardone (1969) was adopted. The Cardone model 
of the marine boundary layer, hereafter referred to as CBL, was an extension 
of the neutral, fixed-terrain model described by Blackadar (1965). The two- 
layer CBL model is schematically illustrated in Figure 2. In the surface 
layer the eddy viscosity (K^^) is a function of height and atmospheric sta- 
bility, according to the Monin-Obukhov similarity theory. Specifying a 
constant equal to its value at the top of the surface layer, and replac- 
ing f by (f + Vg/r), results in a modified Ekman spiral in the upper layer. 
Following Cardone (1969), the surface roughness (Z^) in meters is internally 
determined as a function of the friction velocity 9 



Z . 6.84<10-^) ^ „-3 2 

O U . ^ - 

ic 



4.43(10"^) 



(13) 



The desired boundary layer wind profiles and associated heat flux can 
be determined from the solution of the following set of simultaneous 
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equations , 



Cd =\/2 k sin - a^)/[ln - '|f(L^)] 

2 k Bq sin^ a^/ (()(L*) 



(14) 

(15) 



where 

the geostrophic drag coefficient, 

k = 0.4, von Karman’s constant, 

CX^ = surface cross-isobar ic angle, 

R = V /fZ , surface Rossby number, 
o g o •' ’ 

<t>(L*) = non-dimensional wind shear function, 

= wind profile function, 

= h/L*, stability index, 

h = B^V^/[f + V 0 /r] , height of constant flux layer, 

-4 

B = 3 X 10 , dimensionless constant, 

L' = 0[ln(lO/Z^)- <K10/L' )] /[gk^(T^^-T^^)] , the 

gradient form of the Monin-Obukhov length applied at 10 m, 
ratio of the heat transfer coefficient to the eddy 
viscosity, 

9 = mean potential temperature in the boundary layer, and 

T -T = virtual temperature difference between 10 m and the sea 
va vs ^ 

surface. 

Expressions for and and the iterative method for solving the 

highly implicit set (13) - (15) are given by Cardone (1969). Solutions for 

specify the desired wind profiles and implied heat flux in 

the constant flux layer. It should be noted that the virtual temperature 

difference (T -T ) is used to represent both sensible (Q ) and latent 
va vs s 

(Q ) fluxes from 
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(16) 



H = Qs V L'] 

Individual Q^ and Q^ fluxes are specified according to the ratio (T^-T) 
/(q^-q). The radial wind component was determined by integrating the wind 
profile given by the similarity theory at four points within the constant 
flux layer and by the modified Ekman spiral above* 

Use of the CBL in the atmospheric vortex assumes applicability of the 
similarity theory in wind regimes for which observations are lacking. Thus 
the model results must be critically examined; hopefully, the results will 
stimulate interest in devising and obtaining the measurements necessary to 
improve the boundary layer representation under hurricane conditions. 

Computational aspects . The combined vortex -boundary layer model was 
solved as an initial value problem with centered space and time-differences 
after an initial forward step. Forward time steps were also used periodi- 
cally to reduce the time separation of the solutions. Ocean temperature 
values were constant in space and time. Because the model is to be started 
from non-zero initial wind speeds, temperature and moisture fields that are 
consistent with the vortex intensity must be specified. Although an esti- 
mate of 0 along the radius may be obtained from (2) - (4) using the Vq 
e D 

profile, the values of 9 (or T) and q must also be consistent with the heat/ 
moisture source and the advection. As sufficient data were not available to 
specify all of these variables, two techniques were used to specify initial 
values. In the first experiments the relative humidity was specified and 
the Claus ius-Clapeyron equation and the 0^ definition were solved iteratively 
for the temperature. An alternative approach was to relate the atmospheric 
temperature to the oceanic temperature, and to determine the specific humidity 
from the 0^ value calculated from (2) - (4) . Thus in regions where the 
initial air temperature was too high, the initial q field was too low. If 
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the ocean temperature were fixed, the radial temperature and humidity profiles 
approached the same values within a few hours with both initialization tech- 
niques. Consequently, the second alternative was adopted by specifying an air 
temperature IC less than the ocean temperature. 

RESULTS 

Effect of Ocean Temperature 

A comparison of the maximum tangential wind speed for the two critical 
relative humidity specifications is shown in Table 1 for the two boundary layer 
versions. In both cases there is a systematic increase in maximum Vg as the 
sea surface temperature is increased. Ooyama (1969) also found a change of 
about 20 m sec ^ for a 2C change in sea temperature. When the moisture speci- 
fication dependent on the cloud base is used, minimal hurricane- force winds are 
obtained for an ocean temperature of 28C. With the maximum relative humidity 
specified as 95^, hurricane- force winds are obtained with a sea temperature of 
27c. Similar differences in maximum wind speed are obtained at the same 
ocean temperature for the two moisture specifications, irrespective of the 
boundary layer model. The differences caused by the moisture specification 
decrease slowly with increasing sea surface temperature. 

Riehl (1963) noted the threshold 0 value for producing surface pressures 

e 

capable of sustaining tropical storms and hurricanes was 350K. Parcels with 

surface properties corresponding to 0^ > 350K are required for deep convection 

and warming of the tropospheric column. For higher oceanic temperature, 

the air temperatures and specific humidities throughout the boundary layer 

are increased. Thus the environmental 9 values (at r ) tended to increase 

along with the values near the eyewall. An increase in 9^ at r^ relative to 

0 at r is associated with a pressure decrease as specified in (2), and 
e o 

also with the maximum wind through (4). However, the 0^ of the environmental 



11 



air exceeds 350K for the experiments with high ocean temperatures. If 
only the increase in 9 over the value at r is assumed to contribute to 
the surface pressure drop, the maximum wind speeds shown in parentheses 
in Table 1 result. Only a small increase in wind speed is predicted for 
water temperatures above 30C. Within the framework of three basic assump- 
tions of the model, such an interpretation would suggest a plateau of 
maximum intensity of storms regardless of water temperature in excess of 
a certain value. For a 95^ specification of the critical relative humidity 
at the eyewall, the 60^ values are about 13-15K over the entire range of 
ocean temperatures in Table 1. Wind speeds exceeding the plateau value 
could only be expected if anomalously cool and dry conditions existed 
around a tropical storm over a warm ocean. Such conditions might be ex- 
pected as the storm nears a coastal region, if the inflowing boundary layer 
air originating over land has a sufficiently long trajectory over water 
before entering the central regions of the storm. One can speculate that 
the combination of a favorable environment and the warm coastal waters 
could lead to intensification near the coast, prior to the filling phase 
after the storm center crosses the coast. 

It is somewhat satisfying that the crude representation in terms of 
the bulk aerodynamic transport equations produces maximum intensities which 
are quite similar to those obtained from the more sophisticated CBL. The 
CBL, incorporating the empirical constants as specified by Cardone (1969), 
does appear to be slightly more sensitive to the sea surface temperature 
variations. The similarity of the results is closely related to the criti- 
cal relative humidity specification, since both boundary layer versions 
have large evaporation rates which necessitate the parameterization of con- 
vective clouds. In this sense the external parameters are dominant and 
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other measures must be used to compare the internal parameters of the two 
models • 

Model Parameter Distributions 

One of the objectives of the model formulation was to determine equili- 
brium atmospheric conditions for given ocean temperatures, rather than 
simply specifying fixed air-sea temperature and moisture differences. Values 
of sea-air temperature differences are shown in Figure 3, for the CBL with 
a maximum relative humidity of 95ji. Over most of the region the surface air 
temperature remains essentially constant. A somewhat larger sea-air tempera- 
ture difference is predicted in more intense storms that occur with higher 
ocean temperatures. Near the central region, the air temperature decreases 
rapidly. In the minimum hurricane (27C) the temperature decrease is to a 
value 1C below the water temperature, but for the most intense (31C) storm 
the decrease is more than 4C. Riehl (1963) used a constant difference of 
2.5C for the region with wind speeds exceeding 25 m sec Bell and Tsui 
(1973) found the mean air temperature for 60 ship observations within 
100 n mi of the center was 1.2C less than the ocean temperature. A 5C dif- 
ference was observed as the center of hurricane Hilda passed with 40 n mi 
of the NOMAD buoy (see Marcus and Smith, 1966). Thus, one might conclude 
that the air-sea temperature differences from the model are not inconsistent 
with the variation determined from the few available observations. 

A number of other parameters were examined to evaluate the differences 
between the model results with the two boundary layer treatments. Radial 
variations of several of these parameters are illustrated in Figure 4 for 
two cases in which the wind speed profiles at the top of the boundary layer 
were essentially identical (T^ = 29C in Table 1). 
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In the CBL model the stress is determined from the iterative solution 



for the friction velocity (u^), whereas in the bulk transport equation the 
stress is primarily related to the tangential wind speed at the top of the 
boundary layer. Nevertheless, the surface shearing stress (Tg) profiles 
(Figure 4) were very similar for both models. These values may be compared 
to Tg estimates by Miller (1962) from data collected in hurricane Helene, 
Miller calculated the shearing stress necessary to offset the inward momen- 
tum transport by the radial wind component, which was estimated from air- 
craft data at 800 mb and 200 mb. Average T values were 60, 51 and 28 dynes 
-2 

cm for the radial increments 0-37, 37-74, 74-111 km, respectively. Two 
factors may explain these large T values compared to those from the model. 
Miller assumed the surface tangential wind component was equal to the 800 mb 
value, while a reduction of 0.7 was used in the model. In addition, the 
radial wind component estimates by Miller were larger than the integrated 
v^ values from the model, shown in Figure 4. Thus the transfer to the 
ocean, or the associated radial divergence of momentum, was smaller for the 

model results. Experiments in which the drag coefficient (and thus the 

-3 

heat and moisture exchange coefficients) was constant (C^^ = 2.5 x 10 ) 

did not result in different maximum wind speeds. Even though the heat and 
moisture fluxes at outer radii were increased in these cases, the 0 at the 
eyewall was largely unaffected. Because the maximum wind speed is deter- 
mined by the empirical relations (2) - (4), the primary effect of the larger 
value of was to increase the frictional stress at the ocean surface, and 
this was balanced by increased radial advection of momentum. In this respect, 
the model differs from others (e.g. Rosenthal 1971, Ooyama 1969) in which 
the ultimate intensity reached by the model storm varies inversely with 

Another indicator of the storm intensity is the amount of precipitation. 
In a system such as the hurricane, almost all of the moisture in the inflow 
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layer is precipitated, since the upper level outflow contains little 
moisture. An estimate of the precipitation is then the vertical moisture 
flux through the top of the boundary layer by both the mean and convective 
motions. The moisture source due to both radial advection and evaporation 
is included in this precipitation estimate. An integrated radial component 
of 3 m sec ^ through the 100 km radius with a moisture content of 20 gm kg"^ 
contributes about 12 cm day ^ to the area-averaged precipitation within 
that radius. Since evaporation within 100 km of the center adds only 
another 1 cm day"^ to the precipitation, the total is closely related to 
the strength of the inflow. It might be noted that more than 925^ of the 
vertical moisture flux was by the mean vertical motion. Maximum precipita- 
tion rates, which occurred at the eyewall, for the bulk transport and the 
CBL models compared in Figure 4 were 65 and 30 cm day respectively. 

Heat Extraction from Ocean 

Because the model is eventually to be used to simulate changes in ocean 
thermal structure changes due to the passage of a tropical cyclone, the total 
heat extraction and its radial distribution are important considerations. 
Black and Mallinger (1972) have recently summarized various estimates of 
the heat flux from the ocean. These estimates were based on bulk transport 
equations with surface temperatures, moistures and winds extrapolated from 
aircraft measurements. Maximum fluxes occurred near the radius of maximum 

winds, with the estimates ranging from less than 800 to more than 6000 cal 

-2 -1 

cm day . As Black and Mallinger point out, in many of these flux calcu- 
lations the sea surface temperatures were assumed equal to pre-storm values, 
and thus may overestimate the flux# A second feature of these calculations 
was the rapid decrease in heat flux with increasing radius# Nearly all of 
the flux estimates at five times the radius of maximum winds were less than 
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one-half of the peak values. Similar variations of total heat flux are 
evident in the model results shown in Figure 4, Maximum values occur 
within a few km of the radius of maximum winds, but are considerably smaller 
than estimates based on empirical data. The values in Figure 4 represent 
the total heat removed in the first 24 hours of the calculation, with the 
storm stationary and the ocean temperature fixed. Maximum values of 

-2 

instantaneous latent and sensible heat fluxes at 24 hours were 1.63 x 10 

-2 -2 -1 

and 0.31 x 10 cal cm sec . Thus the Bowen ratio near the eyewall 
was about 20/6. By contrast, near 300 km the Bowen ratio was approximately 
5^. It should be noted that the heat and moisture exchange coefficients 
were equal to 2/3 in the model with the bulk transport equations, since 
it was expected that the molecular exchange of heat and moisture at the 
interface would be inhibited compared to the momentum exchange (Robinson, 
1966). A series of experiments was also run with the heat and moisture 
coefficients equal to the value. Latent and sensible heat fluxes were, 

of course, increased by at least 50^. Extreme values of heat removal 

-2 -1 

ranged from 650 to 6200 cal cm day for ocean temperatures of 27 and 31C, 

respectively. However, the differences in the temperature fields were 

less than 0.5C, with warming at large radii and enhanced cooling near 

the center. The maximum wind speeds were increased by only 3-4 m sec ^ over 

the values shown in Table 1. Associated increases in T^, v and maximum 

U r 

precipitation rate were also noted. It thus appears that the intensity of 
the storm produced by the bulk transport model is rather insensitive to ex- 
change coefficients. Specifying coefficients greater than the minimum 
values necessary to establish the 60^ possible for a particular ocean tem- 
perature produces proportionately larger heat extraction, but not signifi- 
cantly more intense storms. In future experiments with the ocean surface 
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temperatures predicted by the heat balance, the variation of C might be 
expected to have a larger effect. Perhaps one method for determining 
these exchange coefficients will be to determine the oceanic heat loss 
from expendable bathythermograph traces as Black and Mai linger (1972) have 
done for Hurricane Ginger. Calculation of the oceanic heat budget may be 
more reliable than the use of bulk transport equations with asstimed ex- 
change coefficients and with surface winds, temperatures and moistures 
extrapolated from aircraft data. 

The total heat flux after 24 hours in the CBL model was generally 
smaller than that calculated in the bulk transport version (see Figure 4). 
Although accurate observations of the heat flux within the hurricane 
boundary layer are not available, the heat fluxes from the two models 
appear to be underestimates, especially at outer radii. A lower heat 
flux may have been a result of the assumption that the total heat flux 
was proportional to the density (or virtual temperature)gradient near 
the air-water interface. In the mid-latitudes the' stability is primarily 
determined by the large air-sea (or air-land) temperature differences, 
and the moisture content is generally not thought to be important. How- 
ever, in the tropics significant latent heat exchange can occur with 
rather small air-sea temperature differences. Assuming similarity between 
the wind and virtual temperature profiles allows heat transfer even if the 
air and sea temperatures are equal, but it may contribute to an under- 
estimate of the total heat flux. Another factor in the heat flux relation 
[see (16)] is the ratio (a^) of the eddy conductivity (K^j) to the eddy 
viscosity (^) i which is constant if the profiles are similar. However, 
the value of the ratio over the sea is not well known. Thus it has been 
considered a disposable parameter in the model. In a sense this is analogous 
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to the experiments described above in which the heat and moisture coeffi- 
cients were varied relative to Another explanation for the low heat 

flux estimates could be that the u. or stress values are underestimated. 

The key relation for the surface roughness parameter as a function of u^ 
may also be varied to increase the frictional drag at the surface. This 
was the second approach used to test the two- layer model. Point values 
of several parameters are provided in Table 2 from the runs with T^= 29C. 

The surface roughness coefficient (C2 = 0.00428) in (13) and the ratio of 
heat conductivity to eddy viscosity were also varied in runs for the other 
sea surface temperatures, and the resulting changes were similar. 

Enhancement of the heat exchange relative to momentum exchange (Exp.B) 
increased the maximum wind speed about 4 m/sec. This increase in maximum 
wind speed ranged from 3 to 7 m/sec for a T^ range of 27 to 31C. As in 
the case of the bulk aerodynamic model, the increase in surface heat ex- 
change beyond a certain critical value does not dramatically change the 
predicted intensity of the storm for a given T^. The need for observational 
data to determine the proper ratio is again evident. Increasing the 

surface roughness in the model (Exp.C) also resulted in more intense storms, 
with the increase in intensity ranging from 2 to 4 m/sec for T^ values be- 
tween 27 and 31C. As indicated earlier, increases in the surface stress in 
this model tends to increase the supply of energy through radial advection 
without greatly retarding the wind speeds at the top of the boundary layer. 

A general increase in the values of the parameters in Table 2 may be 
expected with the increase in Vg, for example, the larger radial component. 
However, in Exp. C, v^ was increased by about 35^ over the control value, 

compared to a 10^ increase in Exp. B. The stress was about 40^ larger with 

-3 

an implied drag coefficient of 6.4 x 10 near the center in Exp. C. Even 
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at 300 km the drag coefficient was 2.0 x lO”^. The associated inflow 



angles (ranging from 14 to 27 degrees) and ratio of the velocity at the 
top of the surface layer to that at the top of the boundary layer near 
the center were also significantly changed from the control run values. 

By contrast, the above factors related to the predicted wind fields in 
Exp. B were not changed beyond the expected variations due to the in- 
creased maximum wind speed. 

The total heat removed from the ocean in Exps. B and C was larger 
than in the control run. As expected, the increase was about 50^ when 

the ratio was increased from 2 to 3. Peak values ranged from 410 

-2 -1 

and 3940 cal cm day for equal to 27 and 31C, respectively. En- 
hancement of the momentum exchange at the surface produced correspondingly 

larger heat fluxes, with the peak values ranging from 420 to greater than 
2 -1 

4670 cal cm day for the same values. The larger precipitation rate 

near the center in Exp. C is evidently due more to a larger v^ value 
than to the increase in heat gain, because the 50^^ larger heat gain in 
Exp. B is accompanied by only a 20^ increase in precipitation rate. This 
is consistent with the primary moisture source being at larger radii. 

The inflow angles and implied drag coefficients in Exp. C appear excessive, 
particularly at outer radii where more observations are available for com- 
parison. The mechanism through which the heat transfer is enhanced is 
probably in the wind-generated spray (Riehl, 1954), which exposes a much 
greater surface area for heat transfer from the water to the air. It 
appears that the momentum exchange is not correspondingly enhanced by the 
blowing spray (see Monahan 1966, and Wu 1973). If this is the case, the 



roughness. This could be incorporated in the model, and would likely 




function of the amount of spray, or the surface 
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increase the maximum wind speed a small amount. A different representa- 
tion of the heat loss from the ocean versus radius would also result. 

The flexibility of the CBL model, and the fact that future research 
results from air-sea interaction studies will likely be cast in the 
framework of similarity theory, or external parameters such as Z^, is 
thought to be sufficient to justify this more complex and time-consuming 
boundary layer model. 



CONCLUDING REMARKS 

Summary 

An empirical wind profile and an expression for the maximum wind 
speed in terms of the radial increase in 9^ have been used to specify the 
winds at the top of the boundary layer in a S3nmmetrical hurricane. Using 
these winds with the bulk aerodynamic transport equations, and as an input 
to a two-layer, baroclinic marine boundary layer model by Cardone (1969), 
solutions have been obtained for the moisture and potential temperature 
distributions within the boundary layer. The vortex in turn responds to 
the equivalent potential temperature gradient induced by the balance be- 
tween the heat /moisture sources and the advection by the boundary layer 
winds. Equilibrium maximum wind speeds vary from minimal hurricane force 
to 80 m sec ^ for ocean temperatures between 27C and 31C. This might be 
considered the maximum storm for the specified ocean conditions; in actual 
storms the external environmental effects would determine if these wind 
speeds would be achieved. For a given ocean temperature, maximum wind 
speed in the model is shown to be very sensitive to the maximum relative 
humidity of the surface air near the eyewall, which is specified as an 
external parameter in an attempt to parameterize the convective flux through 
the top of the boundary layer. The bulk aerodynamic and the Cardone 
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boundary layer approach give similar maximum wind speeds for the same 
external conditions. There is a certain arbitrariness in the total heat 
flux from the ocean for a given viscosity in both boundary layer versions* 
Observations of heat loss from the ocean are needed to establish the 
appropriate constants in the high wind speed conditions that characterize 
hurricanes. Various parameters predicted by the Cardone model, such as 
the radial mass flow and inflow angles, depth of the boundary layer, verti- 
cal motion and turbulent fluxes, may be evaluated by aircraft equipped 
with inertial guidance systems that are available. 

Possible Applications of the Model 

The original objective was to develop an atmospheric vortex capable 

of simulating the effect of a hurricane on the ocean thermal structure. 

Since the ocean temperature affects the intensity of the vortex, a mutual 

interaction may occur if significant cooling occurs in the ocean. Leipper 

and Volgenau (1970) suggested that the subsurface thermal structure must 

also be considered in forecasting intensity changes in hurricanes. Many 

experiments are planned to evaluate the effect of gradients in both sea 

surface temperature and thermocline depth, A layer ocean model similar to 

that used by Kraus and Turner (1967) is to be tested. 

The adjustment time in the model is relatively short (see Figure 5). 

For T = 31C, the increase in Vq from 30 to 70 m sec ^ occurred within 
w Umax 

11 hrs . Thus the model experiments can be run with a variety of sea tem- 
perature conditions. In the present form the vortex is stationary, so 
that allowing the sea temperature to vary will simulate the rate of decay 
of the vortex due to cooling at inner radii. The method of moving the 
vortex will be based on the integrating effect of strong tangential winds 
in producing properties which are nearly axisymmetric. Application of the 



21 



two-dimensional model along four or more diagonals will be used to esti- 
mate the symmetric component --which determines the intensity of the storm. 
In this way the model vortex may be moved over a varying ocean thermal 
field. One advantage of the present model over previous numerical models 
is the possibility of simulating the symmetric component of actual storms. 
The model can be initialized to a particular maximum wind speed by adjust- 
ing the exchange coefficients or other external parameters, for the 
specified ocean temperature. Experiments are planned to test this capa- 
bility with real data by moving the model storm along the actual track with 
the external parameters held constant. One should then be able to assess 
the predicted changes in maximum intensity of the storm and in the ocean 
thermal structure. 

Finally, the model could be useful in testing proposed hurricane modi- 
fication schemes that are dependent on reduction of the interface fluxes 
of heat and moisture, or the inflow within the boundary layer. It has 
been proposed (Project Stormfury 1968) that spreading of mononuclear films 
could inhibit evaporation from the ocean and thus reduce the energy source 
for hurricanes. Results from the model verify this reasoning, particularly 
if the moisture content reduction could be maintained in the inner region 
of the storm. With currently available organic materials, the film is 
dispersed by the high wind speeds. It might be possible to efficiently 
test many patterns of material spreading to determine the most effective. 
There may be some optimum radius beyond which the materials cannot be 
economically laid down for a given benefit from reducing the maximum wind 
speed. A more recent proposal by Gray (1973) to modify the hurricane by 
carbon dust seeding within the boundary layer might also be tested with 
this model. Gray's hypothesis is that artificial warming of the atmosphere 
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boundary layer due to absorption of insolation by the carbon dust may cause 
an additional mass flux out of the boundary layer. A decrease in the low- 
level inflow is expected to decrease the inner core maximum wind velocities, 
and thus reduce the storm damage. The potential benefit of the various 
hurricane modification schemes is enormous--and so also is the possible 
damage from misguided attempts. Wherever possible, these schemes should 
be tested on a variety of numerical models, including the type proposed here. 
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TABLE !• Maximiun wind speed (in^ec) in the model as a function of sea 

temperature for the two boundary layer versions. The critical 
moisture distribution was by specification of the cloud base 
or by the maximum relative humidity. Values in parentheses 
indicate speed if environmental 0e values exceeding 350K are 
used in Eq. (5) rather than a fixed 6e = 350K at r^. 

Ocean Bulk Aerodynamic CBL 

Temperature Moisture Specification Moisture Specification 

(deg C) Cloud base 75 to 95^ Cloud Base 75 to 95 ^ 



31 


70 


76 


65(51) 


76(62) 


30 


58 


64 


56 (47) 


66(61) 


29 


48 


56 


44 (43) 


56 


28 


36 


47 


32 


46 


27 




36 




33 
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TABLE 2 



. Significant parameters at inner (I), middle (M) and outer (0) 
regions from CBL model upon variation of the surface roughness 
(coefficient C2) and the ratio of heat conductivity to eddy 
viscosity (^^) • Ocean temperature held fixed at 29C. 



Variable 


Region 


A 


Experiment 

B 


C 


C2 


- 


0.00428 


0.00428 


0.025 




- 


2 


3 


2 


VA[m sec ] 

V Im sec J 

"" -2 
Stress [dyne cm ] 


I 


56.5 


60.5 


58 


I 


1.4 


1.5 


1.9 


I 


37 


43 


50 


\ 


I 


.60 


.62 


.50 




0 


• 

00 


• 

00 

00 


.82 


Inflow 


I 


22 


23 


27 


Angle [deg] 


0 


10 


10 


14 


C X 10=' 


I 


2.9 


3.1 


6.4 


u 


0 


1.1 


1.0 


2.0 


Precip [cm day 


I 


33 


38 


51 


Total heat removed 


I 


980 


1460 


1670 


[cal cm ^ day 


M 


180 


100 


115 



27 



LIST OF ILLUSTRATIONS 



Figure 1. 



Figure 2* 
Figure 3. 



Figure 4. 



Figure 5, 



Schematic model of atmospheric vortex and associated inflow 

layer. The exponent x is -1 within r^ and 1/2 between r^ 

and r (see text), 
o 

Schematic of the Cardone (1969) boundary layer model. 

Radial variation of the air minus sea temperatures for the 
CBL with the ocean temperature fixed at 27, 29 and 31C. The 
curves indicate the steady state values with the ocean tem- 
perature held fixed. 

Comparison of radial profiles of selected parameters from 
the bulk aerodynamic boundary layer equations (solid line) 
and the Cardone (CBL) model (dashed). The radial velocity 
at the top of the surface layer (v^), from the CBL model is 
shown in the middle panel. 

Variation of maximum wind speed, starting from an initial 
wind speed of 30 m sec for different ocean temperatures. 
The arrows indicate the time at which the specific humidity 
value at the eyewall first exceeded the critical relative 
humidity. 
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Figure 1. 



Figure 2. 




Schematic model of atmospheric vortex and associated inflow 

layer. The exponent x is -1 within r^ and 1/2 between r^ 

and r (see text), 
o 
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Schematic of the Cardone (1969) boundary layer model. 




Radial variation of the air minus sea temperatures for the 
CBL with the ocean temperature fixed at 27, 29 and 31C. The 
curves indicate the steady state values with the ocean tem- 
perature held fixed. 



Figure 3 



100 
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Figure 4. 



Figure 5. 




Comparison of radial profiles of selected parameters from the 
bulk aerodynamic boundary layer equations (solid line) and 
the Cardone (CBL) model (dashed). The radial velocity at the 
top of the surface layer (v"^), from the CBL model is shown in 
the middle panel. 




Variation of maximum wind speed, starting from an initial 
wind speed of 30 m sec" , for different ocean temperatures. 
The arrows indicate the time at which the specific humidity 
value at the eyewall first exceeded the critical relative 
humidity. 
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